Rayleigh-Plateau instability causes the crown splash 
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The impact of a drop onto a liquid layer and the subsequent splash has important implications for 
diverse physical processes such as air-sea gas transfer, cooling, and combustion. In the crown splash 
parameter regime, the splash pattern is highly regular. We focus on this case as a model for the 
mechanism that leads to secondary droplets, and thus explain the drop size distribution resulting 
from the splash. We show that the mean number of secondary droplets is determined by the most 
unstable wavelength of the Rayleigh-Plateau instability. Variations from this mean are governed by 
the width of the spectrum. Our results for the crown splash will provide the basis for understanding 
more complicated splashes. 



The impact of a drop with a thin film of the same liq- 
uid produces a spray of secondary droplets that results 
from the emission, expansion, and breakup of one or more 
sheet-like jets. Splashes are essential to diverse physical 
processes and applications such as gas transfer across the 
air-sea interface [l[ , cooling Q , coatings [1| , and combus- 
tion 3 . The spatial pattern and size distribution of sec- 
ondary droplets are in general highly complex Q , varies 
qualitatively with experimental conditions 0, 0| , and has 
yet to be understood f§|. Here we focus on the crown 
splash, as exemplified in Edgerton's iconic photograph 
Milk Coronet [9], and show that number of secondary 
droplets is governed by the Rayleigh-Plateau instability. 
Our results can be extended to study more complicated 
splashes such as when there are multiple jets in a single 
impact event or the jets are irregular. 

Figure [T] shows the end stage of a crown splash in 
which the rim of a sheet-like jet breaks into secondary 
droplets distributed almost uniformly along its perime- 
ter. The name of the splash derives from the the resem- 
blance of this final stage to a crown. The events culmi- 
nating in a crown splash begin with a smooth cylindrical 
sheet-like jet shooting outward and upward. The lead- 
ing edge of thisjet is pulled by surface tension towards 
the sheet [lol . , and grows in diameter as it entrains 
fiuid from the sheet. Next, the rim develops a symmetry- 
breaking corrugation, and in a much later nonlinear 
phase of the original instability, the rim's crests sharpen 
into jets which pinch off to form secondary droplets. 

Due to the high speed and small scale structure of 
a splash there are few quantitative time-resolved obser- 
vations, and basic questions regarding the origin and 
evolution of splashes remain unanswered. Leonardo da 
Vinci [I4I appears to be the first to recognize the regu- 
larity of splashes. Worthington was the first to argue for 
a surface-tension-driven instability as the origin of this 
behavior [13|], yet never tested his ideas quantitatively. 
The numerical computations of Rieber & Frohn fl3| sup- 



port a surface-tension-driven mechanism, whereas those 
of FuUana & Zaleski [l^] argue against it. 

Here we present the first experimental demonstra- 
tion that the selection process is governed by a surface- 
tension-driven instability. Our experiments identified the 
parametric regime for crown splashes and measured the 
evolution of these splashes. A 10 cm diameter A/10 glass 
optical flat was placed on the bottom of a 12 x 13 cm 
glass-bottomed container. Fluid was added to the con- 
tainer until the optical fiat was submersed, forming a 
film of height h above the optical flat. The orientation 
of the container was then adjusted so that the fiat lay 
parallel to the fluid's surface to within 3 x lO""* radi- 
ans. The depth of the layer was varied between 150 and 
300 /im depending on the fiuid in order to maintain a 
constant ratio of drop size to layer depth. A single drop 
was released from a gravity fed 30 gauge needle at fixed 
height above the liquid layer at a rate no faster than 
one drop per 10 s, which ensured that the liquid layer 
fully relaxed between splashes. Our data on the morphol- 
ogy of splashes was obtained with a high speed camera 
(Phantom 5.0) viewing the impact from the side. Our 




FIG. 1: Crown splash (silicon oil: Re=966, We=874, H*=0.2). 
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FIG. 2: Qualitative character of impact for H*=0.2. No 
splash (black small circles), crown droplets with (green stars) 
and without (white squares) microdroplets, and microdroplets 
without crown droplets (red diamonds). The solid blue square 
indicates the parameter set for all experiments reported here. 
A crown splash forms exclusively in the Crown droplets regime. 

data for the evolution of the rim was obtained from im- 
ages recorded from directly below through the glass sub- 
strates with an SLR (Canon 20D) and a single 600 ns, 
6 J pulse from a spark flash (Palflash 501, Pulse Photon- 
ics Ltd.). The flash was triggered by the drop cutting 
a laser sheet focused onto a photodiodc. The triggering 
event was reproducible to within ±5 /is. The speed of 
the drop at impact was measured with the high speed 
video camera. The trigger was fed into a delay gener- 
ator which fired the flash. By varying the delay time, 
the evolution of the impact was recreated from a com- 
posite of still images. The virtue of this technique is that 
it produces much higher spatial and temporal resolution 
than can be achieved with a high speed video camera. 
Our measurements of the crown splash were performed 
with silicon oil (density p = 0.92 g/cm'^, surface ten- 
sion 7 = 21 dynes/cm, dynamic viscosity 77=5.2 cP). The 
impact speed and drop diameter were kept constant at 
U = 326 cm/s and D — 0.155 cm, and the layer depth h 
was either 150 or 300 /im. 

Crown splashes occur in a limited range of experimen- 
tal parameters. The dimensionless parameters for de- 
scribing droplet impact in the absence of a surrounding 
gas are the Weber number We = '^^^ , the Reynolds 
number Re = — , the Froude number = and the 

V ' gD ' 

dimensionless substrate fluid depth H* — h/D, where 
g is the acceleration due to gravity. Past studies ig- 
nored the ambient gas, and gravity on the basis that 
the applicable dimensionless numbers are small. We fol- 
low this practice here, and hence only report the We, 
Re, and H*. (The recent work of Xu et al. fli*! found a 
significant influence of the ambient air on drop impact 
on a dry solid. We believe the effect of air to be much 
weaker in our experiment, since there is no moving con- 
tact line.) We measured the phase diagram shown in 



Fig. [2] from a large number of experiments with different 
fluids, droplet diameters, and impact speeds for a fixed 
non-dimensional depth H* — 0.2. Crown splashes appear 
only in the regime labeled Crown Droplets in Fig. [21 Out- 
side of this domain, splashes are more irregular and com- 
plicated 0] . Our measurements of the rim evolution cor- 
respond to dimensionless numbers Re = 894, We = 722, 
and H* = 0.1 ± 0.02 or 0.2 ± 0.01. 

From our images, such as the example in Fig. [31 we 
extracted the average radius of the rim r^, the radial dis- 
tance of the rim from the impact center i?, and the root- 
mean-square amplitude and peak wavelength of the rim 
corrugation. We processed the images to extract the po- 
sition of the the inner and outer edge of the rim to within 
±6 /im. The average radial position for each edge was 
determined from a fit to a circle, and the radius of the rim 
ro{t) and the radial position of the rim R{t) were taken as 
the half the difference and the average of these radii, re- 
spectively. These quantities are plotted in Fig. [H From 
the outer rim data we calculated the power spectrum, 
and the continuous wavelet transform c(0fc , s) , where Ok 
is the translation variable and s is the scale variable, us- 
ing a complex Morlet wavelet. The peak wavelength was 
taken to be the maximum of c(s) = \c{6k, s)p. 

The peak wavelength shifts to larger values for later 
times which corresponds to a coarsening of the corruga- 
tion. Furthermore, the corrugation docs not grow evenly 
at all points on the rim, but rather nucleates at several lo- 
cations. These domains grow and merge, consistent with 
the growth of the unstable modes from random noise. At 
later times the corrugations sharpen and emit droplets. 
The number of droplets is set by the ratio of the rim cir- 
cumference to the peak wavelength, as demonstrated by 
the examples in Fig. [31 and is hence determined by the 
instability. 

We compared our measurements to a theoretical cal- 
culation based on the physical idea that the rim behaves 
like a cylinder of fluid subject to surface tension forces, 
ignoring the effect of the sheet attached to it. Such a fluid 
cylinder is susceptible to the Rayleigh-Plateau instability. 




FIG. 3: Crown splash (Re = 894, We = 722, & H*=0.2) from 
below at t = 1.85 ms and f = 3.15 ms after impact showing 
the one-to-one correspondence between instability wavelength 
and the number of droplets. Yellow arrows define rim radius 
To and rim's radial distance from the impact center R. 
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FIG. 4: Radial distance of the rim from the impact center 
and radius of the rim versus time for Re — 894, We = 722, & 
H*=0.2. 
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FIG. 6: Normalized power spectra versus angular frequency 
for various times after impact for Re=966, We=874, 11=0.2. 
Left: Solid (red) curve calculated from measured rim shaped. 
Dashed (blue) curve calculated from Rayleigh-Plateau model 
using the measured rim radius and circumference. The the- 
oretical peak wavelength is given by the peak of the dashed 
curve. Right: Solid (red) curve was evolved from noisy initial 
conditions using the Rayleigh-Plateau model. Dashed (blue) 
curve is the same as in Left. 



FIG. 5: Ratio of measured to theoretically expected peak 
wavelength versus rim radius for Re=894, We=722, & H*=0.1 
(triangles) or H*=0.2 (squares and circles). Error bars are 
primarily due to the uncertainty in the measured radius of 
the rim. Inset: log-linear plot of the measured root-mean- 
squared amplitude (blue circles), an exponential fit to this 
data (dashed red), and the value computer from the model 
from noisy initial conditions (solid blue) for H* = 0.2. 



which also accounts for the decay of jets mto droplets [17|. 
Its characteristic feature is a periodic variation of the rim 
radius with a period proportional to the average radius 
To- We start at i = with a rim radius r{9j,0) — ro + X 
defined at angular positions 9j , where X is a gaussian dis- 
tributed random variable. Each of the fourier modes of 
the rim grows exponentially with a growth rate am , which 
depends on the reduced wave number qm of the perturba- 
tion, made dimensionless with rg. However, the instanta- 
neous value of qm, and thus am, changes in time accord- 
ing to qm{t) = mro{t)/R{t). Thus the total growth of a 
particular amplitude is determined by the time integral 
over the instantaneous growth rate obtained from clas- 
sical theory. For example, if R{t) increases sufficiently, 
q may become smaller than one and an initially stable 
mode starts to grow. On the other hand, the stretching 
of the rim always decreases the amplitude of the per- 
turbation [13], since longitudinal stretching causes fluid 



elements to contract in the radial direction. This leads 
to a reduction of the amplitude proportional to R{t)^^''^. 
The combined effect gives [l9| : 



r{9,t) = r„ + 
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where qm = "m'-^, P = ^\/iu/ro, cTo = ^/ijrf/p, and 
£i, = v^pI^. The peak wavelength at a any given in- 
stant defined as the mode with the most power, is 
determined by not only the instantaneously most rapidly 
growing mode but also by the history of the other modes. 
Thus the peak wavelength is given by the maximum of 

exp I Jjj* dt' crm(t')|. It bears emphasizing that these cal- 
culations depend on the experimental determined vari- 
ables 7'o(t) and i?(t), and that there are no adjustable 
parameters. 

The ratio of the measured to predicted wavelength is 
plotted versus the rim radius in Fig. [5l The data shows 
that the wavelength is proportional to the expected value 
to within ±20% over a large variation of the rim radius, 
as expected for a Rayleigh-Plateau instability. Our re- 
sults also account for irregularities in the pattern. Under 
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the action of the Rayleigh-Plateau instability, perturba- 
tions away from the most unstable wavelength are am- 
plified as well, albeit at a smaller rate. As shown in 
Fig. [HI the width of the spectrum of a white noise sig- 
nal evolved by our equations is equivalent to that of our 
data. Thus the irregularity is governed by the width of 
the central peak of the spectrum, a phenomenon which to 
our knowledge has never been considered quantitatively 
for the Rayleigh-Plateau instability [3l ■ Instead, typical 
measurements of linear instability impose a wavelength 
from the outside so as to produce a regular pattern, and 
irregularity is often attributed to nonlinear effects [20| . 

The root-mean-square amplitude of the rim deforma- 
tions is poorly reproduced by the our calculation, but for 
understandable reasons. As shown in Fig. [5l[ inset), our 
data is best described by an exponential growth with a 
timescale of 0.59 ms for H*=0.1 and 1.1 ms for H*=0.2. 
The expected root-mean-square amplitude, calculated 
from white noise initial conditions, is not exponential- 
in contrast to a typical linear stability analysis-because 
of the time dependence of the growth rate. On a positive 
note, the observed values of the growth rate fall within 
the range of expected instantaneous rates from the clas- 
sical theory: 0.1-1.2 ms for H*=0.2 and 0.05-0.72 ms for 
H*=0.1. 

We believe the deviations are a result of our simplifying 
assumption that the rim is cylindrical, and the unstable 
mode is varicose. First, even if the sheet is thin com- 
pared to the rim radius, the attachment constrains the 
evolution of the rim and the reduction of surface energy 
due to corrugation is not as great as it would be for an 
unconstrained cylinder. A model calculation in which 
the curve along which the sheet attaches to the rim is 
kept straight, as observed experimentally, shows that the 
growth rate decreases. Second, the gradual transition 
from rim to sheet will slow the growth of the instabil- 
ity. For example, in the extreme case in which the rim is 
simply a semicircular end of the sheet, there no surface 
energy to be gained from a local reduction of the rim 
radius, and thus no instability at all. 

We also considered and rejected other proposed mech- 
anisms to explain the symmetry breaking of the crown. 
Bremond and Villermaux '21] suggested a Rayleigh- 
Taylor instability stemming from the deceleration of the 
rim. Our measurements of the rim's deceleration show 
that the predicted wavelength is 2 times greater at the 
earliest times and 8 greater at late times. Gueyffier & Za- 
leski proposed a Richtmyer-Meshkov type instability due 
to the sudden deceleration during impact [25| . Our mea- 
surements show that the instability grows exponentially 
which is inconsistent with the slower linear growth of a 
Richtmyer-Meshkov instability. Yarin proposed a non- 
linear amplification mechanism [S^l • This mechanism 
does not select a particular wavelength in contradiction 
to our measurements which clearly show wavelength se- 



lection. 

In conclusion, we have demonstrated for the first time 
the selection of a preferred wavelength that sets the num- 
ber of secondary droplets resulting from a crown splash. 
Its number is proportional to the circumference of the 
splash sheet, divided by the radius of the rim. The ir- 
regularity of the crown is set by amplification of random 
noise by the spectrum of growth rates of the Rayleigh- 
Plateau instability. The similarity of the splash patterns 
for other parameter regimes suggests the same mecha- 
nism to be at work, but with a broader range of initial 
perturbations. 

We thank Daniel Bonn for comments on the 
manuscript, and Christophe Josserand for discussions. 
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